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ABSTRACT 



A mathematical model was developed to include line 
effects in the submarine hydraulic system dynamic 
performance analysis. The project was undertaken in an 
effort to demonstrate the necessity of coupling the entire 
hydraulic power network for an accurate analysis of any of 
the subsystems rather than the current practice of treating 
a component loop as an isolated system. It was intended 
that the line model could be coupled with individual elements 
in modular computer programs to facilitate design and alter- 
ation decisions. 

The line model was developed and coupled with a submarine 
driving- system servo-loop in order to demonstrate the above 
premises. However, due to programming difficulties 
additional work is necessary to successfully complete the 
analysis . 
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I. INTRODUCTION 



A. GENERAL DISCUSSION 

The vast majority of hydraulic systems currently in use 
are of a dedicated type. That is, they are designed for a 
specific function with consistent and well defined constraints. 
In most cases they are further limited to a loop in which 
the pump supplies power to a single power element. Such 
systems can be delivered as compact units with a minimum of 
connecting piping. Some examples of dedicated systems in 
Naval use are: surface ship steering units, gun mount 

train and elevation units, and ammunition conveyor systems. 

Although submarines have employed dedicated systems in 
the past, all modern U.S. submarines use constant pressure 
systems [Ref. 1, Tab. 561-1]. The constant pressure systems 
contrast with the dedicated type in that the supply and 
return headers extend the entire length of the ship and serve 
a multitude of power elements. The total piping length of 
these systems may exceed the length of the ship by a factor 
of four or more . 

B. THE NEED FOR A LINE MODEL 

For the packaged configuration of the dedicated system, 
with its relatively short piping lengths , the current practice 
of neglecting line dynamics and lumping the system volume 
has produced acceptable analytical results. However, even 
in systems with low line-volume to total-volume ratios, the 
computed values of the hydraulic natural frequency (u^) 
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are frequently in excess of 407, higher than measured values 

['Ref. 2]. Since sets the upper limit on system response 
and is affected by the volume as well as line dynamics, a 
model which includes these factors would undoubtedly produce 
more accurate results . 

It was the intent of this project to demonstrate that 
the submarine hydraulic system, with its extremely long lines, 
requires that the line dynamics be included to achieve any 
degree of accuracy in analytical solutions. 

C. POTENTIAL FOR SYSTEM MODEL 

Discussions between the author and Mare Island Naval 
Shipyard Hydraulic Design Division personnel (code 260) » 
as well as a Naval Sea Systems Command representative 
[Refs. 2 and 3] , disclosed that there is currently no computer 
program to model the dynamic behavior of a submarine hydraulic 
system. Such a program would be useful in the design of new 
systems, both for sizing of lines and elements and for 
determining the response of components served. For existing 
systems, it would allow a determination of the effects of 
alterations prior to their implementation. 

Additionally, a computer program of a dynamic system 
model would be advantageous in analyzing vibrations trans- 
mitted through the system. The "noise" generated by 
vibrations is a critical element of the "Submarine Sound 
Quieting Program," while the fatigue due to vibrations is 
important to the life of mechanical elements in the system, 
especially the seals [Refs. 5,6, and 7]. 
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The development of a dynamic line model which could be 
programmed in a modular format and coupled with models for 
the individual elements would facilitate the creation of 
such an overall system model. 
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II. SYSTEM MODELED 



A. MAGNITUDE 

It was recognized early in the undertaking that it would 
not be possible to develop a model for an entire submarine 
hydraulic system in the time available. For this reason, 
the stem diving planes servo-loop was chosen for the 
primary development (Fig. 1) . This sub-system was considered 
to be representative of one of the more complex loops in the 
overall system. ( a loop being defined as the fluid path 
through a control valve - to the power element - back 
through the control valve) . Successful modeling of the 
stem planes servo-system would readily permit expansion 
of the system to include the fairwater planes and rudder 
loops in parallel operations. With these three sub-systems 
modeled, the remainder of the hydraulic system could be 
relatively easily developed. 

B. REALISM 

An operational ship's hydraulic system was deemed 
essential as the basis for the model, as it would allow the 
computed parameters to be compared with actual performance 
data. The hydraulic plant of USS PERMIT (SSN 594) was 
chosen as the system to be modeled. It was considered 
representative of a modem submarine hydraulic system and 
the technical information was readily available at Mare 
Island Naval Shipyard. 
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STERN PLANES SYSTEM 



Ill . HYDRAULIC FLUID 



A. TYPE 

The hydraulic oil used in all U.S. submarine internal 
high pressure (3000 psi) systems is Navy symbol 2190TEP. 

This is a petroleum based fluid conforming to military 
specification MIL-L-17331 [Ref. 8] . It exhibits 
approximately the same characteristics as an SAE-30-Detergent 
motor oil. 

B . PROPERTIES 
1. Density 

The density (p) of the oil was obtained from Ref. 9 
in terms of specific gravity (SG) . The values supplied 
were 

SG = 0.896 at 74°F. 

SG = 0.0901 at 60°F. 

The density of water at 60 degrees Fahrenheit is 
-5 3 

9.3475 x 10 (slug/in ) and is assumed to be the reference 
for specific gravity unless otherwise specified [Ref. 10] . 
Specific gravity is defined as the ration of the density of 
the oil at the given temperature to that of water at 60 degrees 
Fahrenheit. 

SG = Po/ Ho (1) 

Solving for the density of oil at temperature (T) 
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yields 



p = 8.4221 xlO~ 3 slug/in 3 at 60°F 

p = 8.3754 xlO" 5 slug/in 3 at 74°F 

2 . Cubical Expansion Coefficient 

The cubical expansion coefficient (a) is the fractional 
change in volume of the fluid due to a change in temperature and 
is used in the determination of density, viscosity (y,v), 
and bulk modulus ($) values. 

Density values arrived at in the previous section 
were used to obtain a value for the cubical expansion coefficient. 

Using the form given in Ref. 11, 

p(P,T) = p o [l + (P-P 0 )/$ - a(T-T o )] (2) 

where the subscript ( Q ) denotes reference values of one 
atmosphere of pressure and 60°F temperature. 

Since the values for the density of the oil were 
obtained at one atmosphere 

P = p o 

equation (2) reduces to 

P - 0 o [1 + a (T - T c )] (2a) 

Inserting appropriate values into (2a) and solving for a 
yields 

a = 3.9638 x 10" 4 °F' 1 

3. Viscosity 

a. Kinematic Viscosity 

Reference 12 specifies that 2190TEP oil shall 
have a kinematic viscosity (v) of 190SSU (Saybolt Seconds 
Universal) at a temperature of 130°F. 
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( 3 ) 



Using the standard conversion [Ref. 14j for 
kinematic viscosity, where the time (t) is greater than 
100 seconds (SSU) 

v = 0.220 t - (135/ t) 

yields 

v = 41.09 cs at 130°F 

v = 92.03 cs at 100°F 

Converting from centistokes (cs) to English units was 

_3 

accomplished with the conversion factor (Ref. 15) of 1.549 xlO 
2 

(in / sec) / cs . Thus the values obtained for kinematic viscosity 
were 

v = 0.06367 in 2 /sec at 130°F 
v = Q.1427 in 2 /sec at 190°F 

b. Absolute Viscosity 

Knowing values of density and kinematic viscosity 

permits the determination of absolute viscosity (p) . 

U - pv (4) 

Using the form presented in Ref. 16, 

V - U o e - A (T - T q ) (5) 

at one atmosphere of pressure, where X is the temperature 

coeffcient of absolute viscosity with units of (°F ^) . 

Combining equations (4) and (5) and then 

solving for lambda gives 

X = In (P 0 v o /pv)/(T-T o ) (6) 

Arbitrarily setting 

T = 100°F and T « 130°F 
o 

in conjunction with values of density, and kinematic 
viscosity at these temperatures yields 

18 



100 



T q = 100 

P o = 8.2886 xlO 



T = 130 F 



-5 



v = 1.427 xlO 
o 



-1 



p = 8.1884 xlO 
>-2 



-5 



v = 6.367 xlO 

Substituting these values into equation (6) produced a 
value of 



X = 2.731 xlO" 2 °F“ 1 

Reference 17 presents the pressure conversion 
for absolute viscosity as 

log 10 (y/y o ) = cP (7) 

However, as the other functions are related by the difference 
in pressure divided by the bulk modulus (AP/B), it would 
seem reasonable to assume that c in equation (7) has a form 
of (c^/ B) • This would then yield an equation for absolute 
viscosity at constant temperature of 

M = y Q e -(c-j^AP/B) (8) 

Lacking any other information, c^ was set equal to one, 
and a final form for absolute viscosity was formulated as 

y = y Q e((P/S) - X(T-T Q )) (9) 

where 



P = P-P 

o 

P = reference pressure (o psig) 

P = system pressure (psig) 

B = system effective bulk modulus (psi) 
T = system temperature (°F) 

T Q = reference temperature (°F) 
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4 . Bulk Modulus 



Although the literature on hydraulic systems indicated 
that the bulk modulus (8) would be supplied with the oil, 
visual inspection of shipping containers and bills of 
lading proved futile. Likewise, the military specifications 
(MIL-L-17331) failed to produce any worthwhile information 
in this regard. A nominal value of bulk modulus for the 
oil was therefore assumed to be 220,000 psi , as recommended 
in Ref. 18. 

C . LIMITATIONS 

In a further effort to maintain the scope of the project 
within reasonable constraints, it was decided to treat the 
temperature as constant. Inferences drawn from the Naval 
Ships.' Technical Manuals previously cited and the author's 
own experience, indicated that a temperature of 130°F was 
a reasonable choice. With this value of temperature, the 
equations for the other properties were simplified. 

From equation (2) 

p = 8.188 xlO" 5 + 8.4221 xl0 _5 (P/6) (10) 

and equation (9) yields 

y » 5.213 xlO -6 e (P/6) . (11) 
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IV . MODEL DEVELOPMENT 



A. SYSTEM OPERATION 

Submarine steering and diving control surfaces (rudder 
and diving planes) are positioned through mechanical linkages 
by hydraulic rams . Each surface has normal as well as 
alternate sources of hydraulic and electrical power to ensure 
reliability. For the purpose of model development, only the 
normal mode of the stern planes was considered. 

The normal mode of operation for the SSN 594 class 
of submarine, employs a position control electro-hydraulic 
servo-loop (Fig. 1). In this type of system (Fig. 2), control 
station operator manually positions a stick to achieve a 
desired angle of inclination of the planes. This produces 
an electrical signal which is summed with the feedback 
signal from the position output to create an error signal. 

The error signal is amplified and transmitted to the torque 
of the servovalve. The torque motor displaces the wand spring 
of the servovalve-pilot, which causes an imbalance in pilot 
pressure, which in turn shifts the main spool allowing the 
servovalve to pass oil to position the ram. As the desired 
position is attained, the error is reduced and the pilot 
returns the spool to a neutral positions, blocking the flow 
of oil. 
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B. SERVO- SYSTEM 



The area of interest for the model development was the 
hydraulic fluid flow in the servo-loop as depicted in Figure 3. 
By inputting the valve spool motion and ram load force as 
explicitly defined functions, the influences of the external 
factors involved in the feedback loops (the servovalve torque- 
motor and pilot valve, and the ram loading) were minimized. 

This permitted the generation of flow and pressure equations 
around the loop with a minimum of complications due to the 
ancillary characteristics. 

1 . Servovalve 

The two servovalves commonly in use in the Submarine 
Force are both manufactured by the Sargent Corporation and 
are designated Sanders Model SV 438-10P and SV 438-15P 
[Ref. 19]. Design and performance of both valves are 
similar. The SV 438-10P which is installed in the (SSN-594) , 
was used for this simulation. Figure 4 is a representative 
drawing of the actual valve employed - a four-land four-port 
element. However, comparison with the valve depicted in 
Figure 3 - a three-land four-port model - will show that 
the flow is essentially the same in both valves. It was 
considered easier to follow the flow paths of the valve in 
Figure 3 . 

All technical data for the servovalve was obtained 
from Ref. 20. 

a. Flow Paths 

The pilot valve flow was disregarded, as it has 
its own supply and return lines. 
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Figure 3 Servovalve -RAM Hydraulic Loop 
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Designating the supply flow as Qj_ n , 
flow as Q r , and the other flows by their respective nodal 
points (Fig. 3), the flow equations were derived. Using 



the form given in Ref. 21, 

£Q ln - JQ out - dv 0 /dt + (V c /B)dP/dt (12) 

the results were : 

^in - Ql + Q 2 (13) 

Q r - Q 3 + Q 4 (14) 

Ql * Q 3 + Q 5 + V^/Sj^ + V x (15) 

Q 2 - Q 4 + q 6 + v 2 p 2 /b 2 + v 2 . (i6) 



Here the volume across the lands of the valve was approximated 
as zero. The superscript (') represents differentiation 
with respect to time. 

Since the valve is manufactured from high 
strength alloys , 




as shown in Appendix A. 

In addition to summing the flows through the 
valve, the flow rate was also defined through the lands for 
three different valve positions. 

i 

(1) Null Positions . In the centered (null) 
position, a critical-center servovalve will have leakage 
past the lands. 1 This flow is described in terms of null 



•For a comprehensive coverage of all terms used in this 
paper, the reader is referred to Ref. 2. 
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flow pressure coefficient (K ) , as 



Q - %k co (ap). 


(17) 


The four node flows for X v = 0 were 




<Jl - % K co< P in- F > 


(18) 


^2 “ % K co< P in- p 2> 


(19) 


<J 3 ’ * K co< P l- p r> 


(20) 


^4 “ * K co< P 2- p r> 


(21) 



with 



See appendix B for values K cq used for the model. 

(2) Valve Spool Off-Centered Left . The valve 
travel (X ) was arbitrarily defined as positive to the left 
(Fig. 3). Employing the orifice flow equation [Ref. 22 ] for 
the fluid passage through nodes one and four, and an empirical 
relationship for nodes two and three, the flow equations were 



determined to be 

Ql - C d A/ ^(Pin-EV/P! (22) 

Q 2 - %V-»(W 2 < p in - p 2> < 23 > 

Q 3 - %K co e-X(X v /X vm ) 2 (P!-P r ) (24) 

Q 4 * C d A/2(p 2 -P r )/P r (25) 

where = 0.61 for full periphery rectangular ports and 

A = TTd x X v (26) 

with d^ defined as the land diameter of the spool valve. 

Since a value of leakage flow was not 



determinable from Ref. 20, the empirical relationships of 
Appendices C and D were developed using Ref. 23 as a guide. 
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(3) Valve Spool Off-Centered Right . The 
equations for the "negative" flow direction were similar to 
those developed in the previous section, with the obvious 
exception that Q 2 and became the principal flows. In 
order to maintain their flows as positive values , was 
squared and included in the radical. 

The form for X less than zero was 



% - %K co e-A(X v /X vm ) 2 (P in -P 1 ) (27) 

Q 2 - C d w / 2x£(P ln -P 2 )/P 2 (28) 

Q 3 = C d «/ 2X 2 (P 1 -P r )/p 2 (29) 

<54 - % K co e -«V X vm> 2 ( p 2- p r> < 30 > 

with 

W = ird 1 (31) 



and all other variables as defined in previous sections. 

2 . Lines 

The line equations were generated with the intention 
of accounting for all factors influencing the pressure and 
flow in the system. The schematic diagram of the actual 
system employed in this development is shown in Figure 5. 

For development purposes , only the odd numbered nodes are 
shown here as the even numbered side equations were identical, 
a. Flow And Pressure 

Disregarding all but the main line piping, the 
flow path from the servovalve to the ram was as depicted in 
Figures 3 and 5 were identified by the same numbers . The 
dynamics of these segments were then defined at each point 
when viewed as in Fig. 6. 
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Figure 5 Stern Planes Hydraulic System Schematic 
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FIGURE 6 LINE NODE IDENTIFICATION 

Point A corresponds to nodal point (5) of Figure 3 
and is considered to be located at the entrance of the pipe 
at its junction with the servovalve. Similarly, point B is 
located at the connection of the two different diameter pipes, 
in the entrance of the downstream pipe. It is equivalent 
to node (7) of Figure 3. Location C is just beyond the end 
of the line segment at the entrance to the cylinder, 
identified as node (9) in Figure 3. 

Utilizing these definitions for the nodal points 
and viewing the direction of flow as always being from A to 
C, the dynamic equations for each element were developed. 

(1) Node Five . The pressure and flow equations 
for node five were derived from Ref. 24. 

q 5 - 0 7 + V 57 f/B 57 + V 57 (32) 

P 5 - P x - P 5 K 5 Q 5 / 2 A 5 (33) 

The double- subscripted variables were employed 
to denote average values through the line segment, while the 
single- subscripted variables represent values at the particular 
node . Thus 
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(34) 



57 



- a 5 l 57 



with L being the actual length of pipe from node (5) 
to node (7) and 

= irD^/4 . (35) 

As can be seen in Figure (6) , the diameter of segment A-B 
(5-7) is constant and equal to the diameter at (5) . The 
total volume and actual length will be developed further 
in later sections of this report. 

The average pressure in the interval was 

P 57 = (P 5 + P 7 )/2 (36) 

with 

P 57 =/P 57 dt. (37) 

The change in volume with respect to time (V) 
was set equal to zero for all line segments , because the 
lines are made from high strength materials which do not 
permit any significant expansion or contraction. The 
exception to this would be in cases where there are flexible 
sections installed ( Aeroquip -type hoses) . There are 
numerous locations in the system where these are installed, 
but not in the servo- loop modeled. In the instances where 
flexible sections are incorporated in the system, an expression 
for the change in volume would have to be developed based on 
the properties of the material used. 

The value of was based on an abrupt change 
in the cross-sectional area of a pipe, and the development 
is shown in Appendix E. As used in the above equations, the 
sign of changes with the direction of flow such that the 
second term of equation (33) will always reflect a pressure 
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loss for the referenced point. Thus, when flow is positive, 
as shown in Figure (6), will be greater than P^, but when 
the flow is reversed (negative in Figure (6)) P^ will be 
greater than P-^. 

Qy is the flow at node (7) and will be defined 
in a subsequent section. P^ is the pressure in the main 
servovalve cavity and is 

Pj_ = J^dt ( 38 ) 

where was introduced in section IV.B.l.a. 

(2) Node Seven . The flow through node seven is 
exactly the same as that through node five, with subscripts 
changed to reflect the appropriate line segment. 

Q? = Q 9 + V79P 79/^79 + V ?9 ( 39 ) 

P 79 = / p 79 dt * 

The pressure at point (7), on the other hand, 
must reflect the losses due to flow in line segment (57) . 

This was accomplished with a combination of the equations 
in Ref. 24 with those in Appendix E: (40) 

P 7 = P 5 - 32u 5 7Le 57 Q 57 /A 5 D5 - l / k* -P 7 K 7 Q 7 / 2 A 5 . 

The second term on the right side of equation 
(40) is the Hagen-Poiseuille term for fully developed laminar 
flow in pipes . The third term accounts for developing laminar 
flow. The final term accounts for the entry /exit loss at 
cross-sectional area changes, as shown in Appendix E. 

All three terms must reflect a loss of 
pressure for flow in either direction. Since Q 57 assumes 
a negative value for flow reversal, it presents no problem. 
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2 

However, Q^y will always be a positive quantity unless it 
is defined as 

Q^7 = Q 57 MQ 57 l, ( 41 ) 

O 

in which case the sign of Q^y will reflect the direction 
of flow. As noted previously, Ky determines the sign of 
the final term. 

2 

Acj is used in the last term of equation 
(40) , becaue the area employed is always that of the smaller 
of the two at a junction. 

Previously undefined variables are 
Q 57 = (Q 5 Qy)/2 (42) 

Le^y = effective length from (5) to (7) . 
Effective line length will be discussed in a later section. 

(3) Node Nine . The pressure at point nine has 
the same form as that at seven 

Pg = Py“32y ygLeygQyg/AyDy - 1 . 1 4 P y g Q y g / Ay - P qKg Qg / 2 A y (43) 

with 

Q79 = Q79 IQ79I ( 4 - 4 ) 

for the reasons discussed in the previous section. 

Flow through this junction is 

Q9 = Q n + Q e x. 

i 

Qn will be defined with respect to ram flow in a following 
segment of the paper. 

Q e i " C el P 9 • (46) 

The coefficient of external leakage (C e ^) is the allowable 
leak rate past the cylinder tail rods to the atmosphere. 
Appendix F derives the values used, based on information 
contained in Ref. 25. 
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There were no volume terms included in 
equation (45) because the transition from the end of the 
line to the cylinder entrance was considered an infinitesimal 
distance . 

b. Secondary Component Effects 

It is readily apparent from a comparison of 
Figure (5) with Figure (3) that the simplified system omits 
numerous valves, fittings, and other miscellaneous elements 
which are installed between the servovalve and ram. Although 
these "secondary " components have a minor influence individually, 
when combined they contribute significantly to both static 
and dynamic characteristics of the system. 

The approach taken in this development was to 
consider them in three categories: branch lines, minor 
components, and major elements. 

(1) Branch Lines . Branch lines are the segments 
of pipe connected to the main line. They are subjected to 
system pressure but do not ordinarily allow fluid to flow. 

Gage lines, hand positioning branches, and capped sections 
were included in Figure (5) . 

These lines contribute only to the total 

system volume. 

(2) Minor Components . The minor components 
were considered to be any elements installed in the system 
which contributed to the effective length and volume of 
the system but were not specifically listed in the list of 
major components [Ref. 26] . 



34 



The elements in this group were ball valves, 
tees, and elbows. Reference 27 was used as the source for 
information on them and Appendix G explains their development 
for use in the overall analysis of the system. 

(3) Major Components. The deceleration valves 
(decel) and the blocking valve were included as major components. 
These were listed in Ref. 26 with specific pressure and flow 
characteristics. The volumes employed were a "best estimate" 
based on diameter and the author's personal experience. The 
effective lengths were as determined from the procedure in 
Appendix H. 

The function of the blocking valve is to 
transfer to an alternate source of hydraulic supply (the 
Vital System), in the event of a loss of main header pressure. 

As pressures in the two systems are ordinarily the same 
(within lOOpsi. of each other) and the valve has internal 
seals, cross-port leakage is negligible. 

The deceleration valve is a cam operated 
full- flow valve which is installed to prevent over- travel 
of the ram. As the ram approaches its maximum limit of 
travel, a collar on the tail rod mechanically trips the 
decel can, shutting the deceleration valve. This permits a 
fluid cushion to arrest the ram travel rather than allowing 
it to a "hard-stop" (metal - to - metal) at maximum velocity. 
Motion in the reverse direction is permitted by flow through 
the internal check valve until the linkage allows the main 
deceleration valve to reopen. 
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(4) Combined Effects. The various elements 



were combined in tabular form for inclusion in the system 
analysis (Tables 1 and 2) . These were input directly into 
the computer program to solve for the various parameters 
as discussed above. 

3. Cylinder and Piston 

The final element in the system is the cylinder and 
piston combination. Here the power is transmitted from the 
fluid to the stem planes by the piston and its associated 
linkages . 

a. Fluid Dynamics 

The pressure and flow equations assumed a form 
similar to those of the lines. 

(1) Flow, The flow in the cylinder was described 
by 

Q 11 = Q if + ^11 + V 11 P 11 /B 11 

where 

^i l = C if^ P ll" P 12^ 

with 

]_ = /P^l^t (^9) 

Qi^ is the internal leakage from the high 
pressure side to the low pressure side of the piston. The 
determination of C^ is shown in Appendix F. Flow was 
defined as positive from node (11) to node (12); thus, when 
P ^2 is greater than P^, will assume a negative value 

as required by the line equations . 
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TABLE I LINE SEGMENT DATA 
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( 50 ) 



V 11 - A U Xr 



is the effective piston area 



— tt(D^ “ 0^)74 



( 31 ) 



where is the land diameter of the cylinder and D r is 

the diameter of the piston tail-rod. 

y 

r is the velocity of the ram. 

V 11 = A 11 L 11 

where 



( 52 ) 





L 11 - L o + X r 


( 53 ) 


and 


L o * ^ L cy - c p> 


( 54 ) 



L c y is the length of the cylinder tp 
is the thickness of the piston. Ram position was defined 



as positive to the left in Figure ( 5 ) , so that motion to the 
left increases while movement Of the piston to the right 
decreases V-^. For this reason, it was necessary to define 



and ultimately 




X 



r 



( 55 ) 



Q 12 “ " Q if " ^12 " V 12^12 /6 12- (56) 

Thus, the various flows in the system would automatically 

assume the proper sign to reflect the direction of flow. 

( 2 ) Pressure. The pressure equation for the 

cylinder and piston arrangement became 

1*11 = ^9 ■^^11^11^11^11^11 “ 1-14 p-j^Q-j^/A^ . ( 57 ) 

There was no entry loss term because it 
was taken into account with P^. The other difference, with 
regard to the previous line loss equations, is that D^is the 
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land diameter of the cylinder while is the effective 
area of the piston. The effective area was used because 
Q/A is a velocity term employed in defining viscous loss. 

On the other hand, was used to account for the surface 

area upon which the fluid is '"dragging"; i.e. , the walls of 
the cylinder. Since the tail rod is traveling at essentially 
the same velocity as the fluid, it contributes no losses 
due to viscosity. 

As in previous sections, it was necessary 



to define 



Qll " Qu-IQ U I (58) 

to account for flow direction, 
b. Piston Dynamics 

Ram motion was defined in terms of position (X ) , 
velocity (X r ) , and acceleration (X r ) . 



X r - 


/x r dt 


(59) 


• 

X = 
r 


Jx r dt 


(60) 


X r * 


( A ll P ll' A 12 p 12' m oi A'V /M T 


(61) 



The load force (F^) , as described at the outset, 
was input as a ramp function dependent on ram position and 
velocity. 



The flow of oil (m Q ^) term was included to 
account for the force required to maintain. the momentum of 
the oil. It was determined to be 

™oi]_ = (3 12^12 f° r •^r > ^ (62) 



and 



m oil _P 11 V 11 



for X r <0. 



(63) 
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The total mass of the system (Mp) was defined 

as the mass of the ram plus the mass of the oil from the 

low pressure side of the piston to the return tank, such that 

tank 

Hr. - m + l ppVp (64) 

ram 

where the subscript (p) indicates each main line segment 
(flow path) from the ram face to the return tank. 

C . PROGRAM 

The foregoing equations were re-arranged to yield the 
form shown in the attached computer program. They were 
modeled in a Continuous System Modeling Program (CSMP III) 
in the belief that it was the best modeling method available. 
However, numerous attempts failed to accomplish the desired 
result . 

The algebraic equations form a loop for which no solution 
was found. 
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V. SUMMARY 



A. CONCLUSIONS 

It is felt that the development of the line equations 
demonstrates the necessity of including these factors in 
any in-depth analysis of the submarine hydraulic system. 

It is also the author's belief that the approach taken in this 
presentation will yield worthwhile results for any system 
employing extensive quantities of piping. 

As stated in the introduction, a generalized computer 
program of this type would have many benefical aspects for 
both the submarine designer as well as the operator. 

B . RECOMMENDATIONS 

The first recommendation is a possible means of clearing 
the algebraic loop problem currently preventing the program 
from executing. This would require clearing all of the 
pressure and flow equations of reference to any variables 
except for themselves or ones which directly invoke an 
integral. That is 

Py = f (Pyg > P5 y > P^ » Py > £ tC • ) . 

However, it could not take the form where Py was set equal 
to Pg or any of the flow terms (Q^), as they do not invoke 
an integral function. 

The other possibility would be to locate another program 
which could efficiently handle both integral as well as 
simultaneous algebraic equations. 



42 



APPENDIX A 



Determination of Servovalve Volume 



No values were obtainable for the volume of the 
servovalve. Therefore, the following method was employed 
to obtain a realistic approximation. 

Using nominal dimensions given in the technical manual 
[Ref. 20 ] and the knowledge that the valve is made of 
aluminum alloy, a volume to weight ratio was determined. 
Dimensions 18.5 x 8.0625 x 9.5 in. 

Specific Gravity of alloy approximately 2.8 
Given weight 69 lb. 

Density of water at 60°F. 62.41b/ft 3 

V body - 1417 in3 

Y alloy = 0.101 lb/in 3 

W = Vy = 143.3 lb for a solid block 



then 



W given /w = °' 48 



Thus the valve body weighs approximately half of what a 

solid block would weigh. It was therefore assumed that 

3 

the volume of metal in the valve was 708in . Further, 

estimating that the passages for the pilot system occupy 

half of the void area yielded a main valve volume of about 
3 

350 in , and since this volume is that occupied by both the 
supply and the return segments of the valve, half was assigned 
to each. 

V 1 = V 2 = 175 in 3 
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APPENDIX B 



Null Flow Pressure Coefficient Determination 

The allowable null flow leakage was given as 600 cc/tnin. 

at 1200 psi. [Ref. 20]. This was converted to standard form 
3 

of (in /sec) /psi. as follows 

K cq = 600(cc.min) /1200 psi x 1 in 3 /16.39cc x 1 min/60 sec • 
This yielded 

K cq = 5.0844 xl0~ 4 (in'Vsec) /psi. 

and 

%K cq = 2.5422 xlO -4 (in 3 /sec) /psi . 

APPENDIX C 

Empirical Valve Leakage Coefficient 
The leakage flow across a valve port decreases as 
the valve moves to block the port. Since a plot of leakage 
flow (Q^) versus valve position (X v ) is symmetrical, with an 
exponential decay as shown in Fig. C-l, and exact data was 
not given for leakage other than the null position an 
empirical expression was developed. 




Figure C-l Leakage Flow vs Valve Position 
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The most appropriate expression to fit this curve took 



the form of 

V<W e-^V X ™ax ) 2 

where 



= leakage flow at any valve position 



Q^ max = maximum leak rate = %K cq 
X vmax = max i- rau m valve travel from null. 
Representative values were then chosen such that 



V<W - - 02 ac V x ™a x - ' 25 

That is, at 25% of maximum valve, travel leakage through 
the "blocked" port of the valve decreases to 2% of the null 
position leakage. 

Solving for X 

x - tww'VWfVW 2 

X = 62.6. 



In reality X could vary from near zero to about 100. 
Realistically it would have a value between 50 and 100. 

X = 65 was chosen as the value to employ in this 
presentation. 

The form used was 

q £ = %K co e- X <V X vmax > 2 

with the determination of maximum valve travel (X ) shown 

vm 

in Appendix D. 
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APPENDIX D 



De Germination of Maximum Valve Travel 



Reference shows that 

Qmax = 80 GPM = 80 x 231/60 in 3 / sec 
at 



T = 115 F 



AP total = 1200 psi 

Ap port ' 600 psi 
as depicted in Figure D-l. 




6oo psi 



Figure D-l Valve Test Flow Path 



Using the orifice equation 
Q = C d A(2Ap/p)^ 

with 

C d = 0.61 
p = 8.442 xlO" 5 
A =7r ^]_^ v where d^ = 1 in. 

X v - X vm 
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such that 



X vm ’ <V" C d d l< 2 A ye>« 

resulting in X = 0.0426 in. 

This value was used in spite of the fact that it seemed 
large for a “standard 0 pilot controlled servovalve. 



APPENDIX E 



Loss Coefficient Due To Abrupt Area Changes 

The standard form of the pressure loss due to abrupt 
piping diameter change is given as 
H 1 = KV 2 /2g 

where 

is head loss, expressed in feet of fluid 

Y is fluid velocity 

g is acceleration due to gravity • 

Head loss can also be expressed as 
H 1 = Ap/ a 

Y = Pg 

and 

Y is weight density 
while velocity is 

V = Q/A = 4Q/ttD 2 
combining the above equations 
AP = 8 P KQ 2 /tt 2 D 4 . 

K is defined from Figure D-l 
K = %(1-D 1 /D 2 ) for +Q 
K = -(1-(D 1 /D 2 ) 2 ) 2 for -Q. 

NOTE D-^ is always the smaller of the two diameters being 
considered. Q+ 



Figure E-l Flow at an Abrupt Area Change 
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APPENDIX F 



Seal Leakage Coefficients 



Seal leakage was based on permis sable values listed in 
Naval Ship's Technical Manual, [Ref . 25.] The values used in 
this presentation were for '’'used''' seals with maximum allowable 
leakage before replacement is required. This was felt to 
be reasonable since the operating temperature of the oil 
is normally greater than the test temperature, which will 
produce slightly higher lead rates. 

The leak rates were given as one milliliter per ten 
cycles per inch of rod diameter for external leakage, and 
five milliliter per five minutes per inch of rod diameter 
for internal leakage. These were based on normal operating 
pressures and a fluid temperature of 100°F. 

Operating pressures were assumed to be 
P e ^ = 2000 psi 
P^ = 1000 psi 

where P £ ^ is the pressure for external leakage and P^ is 
the pressure for internal leakage. A cycle was defined 
as the time to go from neutral (zero plane angle) to full 
rise to full dive and back to neutral position. For the 
system modeled, stem planes travel is plus or minus 27°. 

Thus, a full cycle would be four times this value, or 108°. 

The normal rate of travel is about five degrees per second. 
Thus, 108/5 = 21.6 sec/cycle or approximately 
3 cycles/min. 
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Since the tail rod used, had a diameter of four inches, 
the allowable external leakage was four milliliters per ten 
cycles . 

Then 



4ml x lcc x 3_cx . i.zcc/ndn 
10 cy 1ml 1 min 



and 

3 

C s = 1.2 cc v 1 min Y 1 in v 1 

el mn * KTZZc 1'673'9'c'c IMtfpiT 

C e £ = 6.1 xlO”^ (in^/sec) /psi 

Similarly the internal leakage for a land diameter of 
8.25 inches, rounded to nine inches, was nine cubic 
centimeters per minute . 

And 

c - 9 

C U 60 x 16.39 x 1000 

C i£ = ^ *15 xlO ^ (in^/sec) /psi 
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APPENDIX G 



Minor Component Analysis 

The values used in the derivation of the minor component 
effects, considered herein, were obtained from Crane Technical 
Paper No. 410, [ Ref. 27] . 

The volume was determined to be 
V = AL 

where A was the cross-sectional area of the pipe entering 
or leaving the element and L was a value determined for 
each type of element. The values for L were chosen as 
Ball valves L = 3 x D 

Tees L = 2 x D 

Elbows L = 1.6 x D 

In cases where the element has differing entry and exit 

diameters, D is the larger of the two. 

The values for the effective lengths were 

Ball valves L = 3 x D 

e 

Flow through tees L g = 20 x D 

Branch flow tees L = 60 x D 

e 

45° elbows L = 16 x D 

e 

90° elbows L = 30 x D 

e 

where the diameters were as above. An additional factor of 
K (appendix E) would have to be included for diameter 
differences in the main flow path. 

The value of effective length cited above, were determined 
for Reynolds numbers (R g ) greater than 1000. Therefore, 
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it was necessary to modify them 



Since 

and 

with 

then 

If 

so that 
or 

and for 



R g = 4Q/ (ttDv) 

L ,= L R /1000 for R <1000 
ed ec e e 

the desired effective length 
L ec = the calculated effective length from above 

L e = L ec (^Q/(ttDv))x10' 3 

R > 1000 then L = L 
e - e ec 

for 4 Q/ttDv>1000 

Qp/Dy> 785.4 L = L 

x - e ec 

Qp/Dy<785 . 4 L e = 1.273 xl0 _3 L ec Qp/Dp 
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APPENDIX H 



Effective Length of Major Valves 



The standard head loss equation 
o 

= K v /2g was rearranged in the form 
AP/pg = KQ 2 /2gA 2 
This yielded 

K = 2aPA 2 / p Q 2 

since 

A =ttD 2 /4 
K =tt 2 APD 4 /8 P Q 2 
and since 

K = L e /D 
L e =tt 2 aPD 5 /8 p Q 2 

To allow Q to be entered in the form presented (GPM) , 
a conversion was included. 

231 in 3 /gal x 1 min/ 60 sec 
so that 

Q (GPM) x 231/60 = Q (in 3 /sec) 

and L e = .0832 PD 5 /pQ 2 
where 

D = diameter of valve (in) 

P = given pressure loss (psi) 

Q = given flow (GPM) 

P = density of fluid at given temperature (slug/in ) 
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APPENDIX I 



TO EMPLOY: 

SUPPLY LINES; D1 IS ALWAYS THE POINT CLOSEST TO THE 
PUMP - THE UPSTREAM POINT 

RETURN LINES; Dl IS ALWAYS THE POINT FARTHEST FROM 
THE TANK - THE UPSTREAM POINT 

I 

DUAL DIRECTION (REVERSING) LINES; Dl IS THE POINT 

NEAREST THE CONTROL ELEMENT (VALVE), VIEWED AS IF 
THE LINE ALWAYS SUPPLIES PRESSURE TO THE POWER 
ELEMENT 

ENDMAC 

INITIAL 

CONSTANT KCCCC=600. ,FT=1200. ,LAMB=65. 

VALVE LEAK RATE OF MAIN SPOOL 

NOTE: LAMB IS STRICTLY EMPIRICAL AND BASED ON AN 
AVERAGE CURVE OF LEAK-RATE VS VALVE DISPLACEMENT. 

CONSTANT XVM=0 . 0425, DLV=1 . 0 , CD=0 . 6 l , Vl=175. , V2= 175. 

CONSTANT 0D5=1.90,TH5=0.25,L5=48.0,NTB5=0. ,ntt5=o., ... 

NL45=0 . , NL95=1 . , NV5=0 . 

CONSTANT OD 6 =l. 9, TH6=0. 25, L6=48. 0, NTB6=0. ,NTT6=0. , . . . 

NL46=0 . , NL96=1 . , NV6=0 . 

CONSTANT CD7=2 . 0 , TH7=0 . 27 , L7=228 . 0 , NTB7=2 . , NTT 7=1 . , NL47=0 . , . . . 
NL97=6. ,NV7=1. 

CONSTANT 0D8=2 . 0 , TH8=0 . 27 , L8=216 . 0 , NTB8=3 . , NTT8=0 . , NL48=0 . , . . . 
NL98=5.,NV8=1. 

CONSTANT 0D71=0 . 75, TH71=0 . 109 , L71=8 . 0 , NTB71=1 . , NTT71=C . , . . . 

NLl+7 1=0 ., NL9 71=0 ., NV7 1=0 . 

CONSTANT 0D72=0. 375, TH72=0. 058, L7 2=120. 0,NTB7 2=0. ,NTT72=0. , . . . 
NL472=0 . , NL972=1 . , NV72=1 . 

CONSTANT 0D73=1.5,TH73=0.203,L73=12.0NTB73=0. ,NTT73=0. , . . . 

NL47 3=0 . , NL97 3=0 . , NV7 3=0 . 

CONSTANT 0D8l=0.75,TH8l=0.109,L8l=8.0,NTB8l=l. ,NTT8l=0. ,. . . 
NL48l=0. ,NL98l=0. ,NV8l=0. 

CONSTANT 0D82=0. 3750, TH82=0. 058, L82=120. 0, NTB82=0. ,NTT82=0. ,.. . 
NL482. ,NL982=2. ,NV82=1. 

CONSTANT 0D83=1.5,TH83=0.203,L83=12.0,NTB83=0. ,NTT83=0. , . . . 

NL48 3=0 . , NL98 3=0 . , nv8 3=0 . 

CONSTANT PDR7=60 . 0 , QDR7=50 . 0 , RDR7=8 . 68E-05, VDV7=60 . , DDV7=1 .094 
CONSTANT PBR7=50.,QBR7=50. ,R3R7=8.68E-05,VBV7=20. ,DBV7=l-46 
CONSTANT E=3.0E07,P0IS=0.3 
CONSTANT B0=200000 . 

CONSTANT PR=50. ,PIN=3000. ,FA=0.1,0DRE=2.0,THRE=0.l48,VR-2750. 
CONSTANT 0ILCC=9.0,CELCC=1.2,PCEL=2000. ,PCIL=1000. 

CONSTANT DLR= 8 . 255 , DRR=4»0 , 0DR=11 . 5, LCY=39 . 875, TR=5. 125 , MR=12 . 1 
CONSTANT PI=3. 141593 
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KCC=0 . 5*KC0 CC/ ( PT*98 3.4) 

THE VALUE OF (0.5) IS USED BECAUSE THE TECH. MANUAL 
LISTS A LEAK- RATE FOR THE WHOLE VALVE; THE 
VALUE SOUGHT IS THAT FOR HALF THE VALVE (ONE SET 
OF PORTS). 

KCOCC IS TECH. MANUAL LEAK-RATE IN (CC/MIN). 

PT IS TEST PRESSURE USED TO ESTABLISH THE LEAK-RATE (PSI). 
KCO IS LEAK-RATE IN ( (lN*-*3/SEC)/PSI) . 

AG=CD*PI*DLV*SQRT ( 2 . ) 

AG IS THE PRODUCT OF THE CONSTANT COEFFICIENTS USED IN 
IN THE SERVO- VALVE EQUATION. 

CFA=FA/1.4 

FA IS THE FRACTION OF AIR ENTRAPPED IN THE- SYSTEM, 

EXPRESSED AS A PERCENTAGE OF THE TOTAL SYSTEM VOLUME. 
CFA IS THE COEFFICIENT OF THE AIR TERM IN BULK MODULUS 
CALCULATIONS. 

k5p , k5m=ent co ( di , d5 ) 

K6P , K6M=ENT CO ( D2 , D6 ) 

K? P , K7M=ENT CO ( D5 , D7 ) 

K8 P , KSM=ENT CO ( D6 , d8 ) 

K9P , K9M=ENTC0 ( D7 , Dll ) 

K10P , K10M=ENT C0(D8,D12) 

THESE COEFFICIENTS ACCOUNT FOR EXIT LOSSES ONLY: 

ENTRY LOSSES ARE CONSIDERED AS EXIT LOSSES FOR 
THE NEXT ELEMENT UPSTREAM. 

vp5,a5,lep5,d5la5,bcs5=vale(pi,e,pois,od5,th5l5,ntt5,ntb5, . . . 

nl45,nl95,nv5) 

VP6 , A6 , LEP6 , D6LA6 , BCS6=VALE ( PI , E, POIS , 0D6 , TH6 , L6 , NTT6 , NTB6 , . . 

NL46,NL96,NV6) 

VP7 , A7 , LEP7 , D7 , LA7 , BCS7=VALE ( PI , E, PO IS , 0D7 , TH7 , L7 , NTT7 , NTB7 » • 

NL47,NL97,NV7) 

vp8,a8,lep8,d8,la8,b(s8=vale(pi,e,pois,od8,th8,l8,ntt8,ntb8, . 

NL)+8,NL98,NV8l) 

V&B1,A7B1,LE7B1,D7B1,LA7B1,BCS71=VALE(PI,E,P0IS,0D71,TH71, . . . 

L71,NTT71,NTB71»NL471,NL971,NV71) 
V7B2,A7B2,LE7B2,D7B2,LA7B2,BCS72=VALE(PI,E,P0IS,0D72,TH72, . . . 

L72 , NTT72 , NTB72 , NL47 2 , NL972 , NV72 ) 
V7B3,A7B3,LE7B3,D7B3,LA7B3,BCS73=VALE9PI,S,POIS,OD73,TH73, . . . 

L7 3 , NTT7 3 , NTB7 3 , NL47 3 , NL9 7 3 , NV7 3 ) 
v8bi,a8bi,le8bi,d8bi,la8bi,bcs8i=vale(pi,e,pois,od8i,th8i, . . . 

L 8 1 , NT T8 1 , NTB8 1 , NL48 1 , KL98 1 , NV8 1 ) 
V8B2,A8B2,LE8B2,D8B2,LA8B2,BCS82=VALE9PI,E,P0IS,0D82,TH82, . . . 

L82 , NTT82, NTB82,NL482, NL982, NV82) 



v8B3,A8B3,LE8B3,D8B3,LA8B3,BCS83=VALE(PI.POIS,OD83,TH83,... 

L8 3 , NTT8 3 » NTB8 3 , NLl+8 3 , NL98 3 , N V8 3 ) 



B0I=1.0.B0 

BCII=2. 0»(1. 0+POIS)/E 

BI1=B0I+BC1I 

BI2=BI1 

B15=BI1 

BI6=BI1 

BI7=BI1 

Bl8=BIl 

BI9=BI1 

BI10=BI1 

BI11=BI1 ■ 

BI12=BI1 

BCDV8=BI1*VDV8 

BCDV7=BI1*VDV7 

bcbv8=bii*vbv8 

BCBV7=BI1*VBV7 

bi57=boi+bcs5/v57 

BI68=B0I+BCS6/v68 

Bl79=BOI+(BCS7+BCS71+BCS72+BCS73+BCS7i4-+BCDV7+BCBV7)/V79 
Bl8l=B0I+ ( BCS8+BCS81+BCS82+BCS8 3 +BCS 8 I 4 .+BCDV 8 +BCBV 6 ) /V*l 

cq57 » cq57s=linco ( a5 , d5 , 0 d5 ) 

CQ58 , cq68s=linco ( A6 , D6 , 0D6 ) 

CQ79 , CQ79S=LINCO ( A7 , D7 , 0D7 ) 
cq8i , cq8is=linco ( a8 , d8 , od8 ) 

CQ11 , CQ11S=LINC0 ( All » Dll » ODR ) 

CQ12 , CQ12S=LINC0 ( A12 , D12 , ODR ) 

PlZ=(PIN-PR)/2,0 

P2Z=01Z 

P57Z=P1Z 

p68z=PlZ 

P79Z=P1Z 

P81Z=P1Z 

P11Z=P1Z 

P12Z=P1Z 



Q7Z=0.0 

q8z=o.o 

Q9Z=0 .0 

Q10Z=0.0 

Q11Z=0.0 

Q12Z=0 .0 

Q57Z=0.0 

Q68Z=0.0 

Q79Z=0.0 

Q8lZ=0.0 

A11=PI*(D11->D11-D11R-x-D11R)A.O 
A12=PI* (D12*D12-D12R*D12R) A.O 
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D11=DLR 

D12=DLR 

D11R=DRR 

D12R=DRR 

A5S=a5-::-a5 

A6S=a6*a6 

A7S=A7*A7 

a8s=a8*a8 



L57=LA5 

a68=la6 

L79=LA7 

L8l=LA8 

lep57=lep5 

LEP7 9=LEP7+LDV7 +LBV7 

LEP68=LEP6 

LEP8l=LEP8+LDV8+LBV8 

V57=VP 5 
V68=VP6 

V79=VP7+V7B1+V7B2+V7B3+VDV7+VBV7 

V81=VP8+V8B1+V8B2+V8B3+VDV8+VBV8 

LDV7=SPECV ( PDR7 , QDR7 , DDV7 ) 
LDV8=LDV7 

LBV7=SPECV ( PBR7 , RER7 , DBV7 ) 

LBV8=LBV7 

VBV8=VBV7 

VDV8=VDV7 



DRE=0DRE-2.0*THRE 

BIR=B0I+2.0*( (l.+POIS)*ODRE*ODRE+(l.-POIS)*DRE*DRE)/. . . 
( E-* ( ODRE-”-ODRE- ORE*ORE ) ) 

CEL=CELCC /(PCEL*983.1^) 

CIL=CILCC/ ( P CIL-ii-98 3 • ) 

XRZ=0 .0 
ZRDZ=0.0 

Dl=3«'”*DLV 

D2=D1 



LC=0 . 5«- ( LCY-TR ) 



DYNAMIC 
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** PROCEDURE SERVO- VALVE *** 



PROCED Ql,Q2,Q3,Q4=SERVO(AG,KCO,LAKB,PIN,PR,Pl,P2, ... 

XV , XVM , RHOl , RH02 , RHC R ) 

CEX=KCO*EXP ( -LAMB* ( XV/XVM)*-»2 ) 

IF ( XV ) 10 , 20 , 30 
10 Q1=CEX*( PIN-PI) 

Q2=AG-::-SQRT ( XV*XV* ( P I N- P2 ) /RHO 2 ) 

Q3=AG*SQRT ( ZV*ZV* ( Pl=PR ) /RHOR ) 

Q4=CEX*(P2-PR) ' 

GO TC i+0 

20 Q1=CEX*(PIN-P1) 

Q2- C EX* (PIN-P2) 

Q3=CEX*(Pl-PR) 

Q4=CEX*(P2-PR) 

GO TO 1+.0 

30 Ql— AG*XV*SQRT ( ( PIN-PI ) /RHOl ) 

Q2= CEX* ( P I N- P 2 ) 

Q3=CEX*(Pl-PR) 

Q4=AG(XV(SQRT( (P2- PR) /RHOR) 

1|0 CONTINUE 
ENDPRO 

** PROCEDURE KPRO *** 

PRO CEED K5 , K 6 , K7 , K 8 , K9 , KlO=KPRO ( Pi , P2 , P$, P 6 , P7 , P 8 , P9 , P10 , K5P , . . 

K5M,K6P,K6M,K7P,K7M,K8P,K8M,K9P,K9M,K10P,K10M) 





IF(P1.GE.P5)G0 
k5=k5m 
GO TO 101 


TO 100 


100 


k5=k5p 




101 


IF(P2.GE.P6)G0 
K 6 =k 6 m 
GO TO 103 


TO 102 


102 


K6=K6p 




103 


IF (p5.GE.P7)GO 

K7=K7M 

GO TO 105 


- TO 104 


104 


K7=K7P 




105 


IF(P 6 .GE.P 8 )GO 

k=k8m 
GO TO 107 


TO 106 


106 


k8=k8p 




107 


IF(P7.GE.P9)GO 
K9=K9M 
GO TO 109 


TO 108 


108 


K9=K9P 




109 


IF(P8.GE.P10)G0 
K10=K10M 
GO TO 111. 


TO 110 


110 


K10=K10P 




111 


CONTINUE 




ENDPRO 
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*** PROCEDURE EFFECTIVE LINE LENGTH 



PROCED LE57=PR057 ( LEP57 , L57 , Q57 , D5 , RH057 , MU57 ) 

CLE57=q57*RH057/ ( mu57*d5 ) 

IF( CLE57.LT. 785. i;)GO TO 200 

le57=l57+lep57 ' 

GO TO 201 

200 LE57=L57+1 . 27 3E-03*CLE57*LEP57 

201 CONTINUE 
ENDPRO 

PRO CED LE68=PR068 ( LEP 68 , L 68 , Q 68 , D6RH068 , MU 68 ) 
CLE68=Q68RH068/ (MU68 *d6 ) 

IF (CLE68.LT. 785. 4) GO TO 202 

LE68=L68+LEP68 

GO TO 203 

202 LE68=L68+1.273E-03*CLE68*LEP68 

203 CONTINUE 
ENDPRO 

PROCED LE79=PR079(LEP79,L79,Q79,D7,RH079,MU79) 

CLE7 9=Q7 9-::-RH07 9/ ( MU7 9-*D7 ) 

IF(CLE79.LT.785.i;)GO TO 201; 

• LE79=L79+LEP79 

GO TO 205 

201; LE79=L79+1. 273E-03-::*CLE79-::-LEP79 

205 CONTINUE 
ENDPRO 

PROCED LE8l=PR08l ( LEP 8 l , l 81 , Q8l , D8 , RH081 , MU8l ) 
CLE 8 l=o 8 l*RH 08 l/ (MU 8 l*D 8 ) 

if(cle8i.lt785.1;)go to 206 

LE8l=L8l+LEP8l 

GO TO 207 

206 LE8l=L8l+l . 27 3E-03-* CLE8l-x-LEP8l 

207 CONTINUE 
ENDPRO 

***. PROCEDURE MASS Op OIL 

PROCED MDOIL ,MT=MOPRO ( MR , XRD , RHOll , RH012 , Vll , V12 , VllD, 
V12D, RK057 , RHO 68 , RH079 , RH08l , RH01 , RH02 , VI , V2 , V57 , V63 , 
V79,V81,... 

RHOR, VR) 

IF(XRD.GE.O.)GO TO i;00 

K0IIrsRH0ll»Vll+RH079*V79+RH057*v57+RH0*Vl+RH0R*VR 
GO TO 1;01 

l;00 M0IL=RH0l2*Vl2+RH08ls-v8l+RH068-:s-RH02*::-V2+RH0R---VR 

MD0IL=RH0l2*V12D 
1;01 MT=MR+MOIL 

ENDPRO 
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PlD=(Ql-Q3-05)/Vl*BlI) 

P2D=(Q2-Q1j.-QS)/(V2*B2I) 

p57D=(qS-Q7)/(b57I*v57) 

P68d= ( Q6-Q8 ) / ( B68l*v68 ) 

P79D=(Q7-Q9)/(B79I*V79) 

P 8 1 D= ( Q8 - Q10 ) / ( B8 1 1-::- V8 1 ) 

PllD=(Qll-OlL-VllD)/(VllI-rrVll) 

P12D=(Q12+QIL+V12D)/(B12I*V12) 

P5=2.0*P57-P7 

P7=2.0*P79-P9 

P9=(Q9-QH)/CEL 

p6=2.0-*p68-p8 

p8=2 . 0-"-p8l- PlO 

P10=(Q10-Q12)/CEL 

s- LINE AND RAM FLOWS *** 

q5=a5*sqrt ( 2 . o* ( pi - p5 ) / ( rho5*k5 ) ) 

q6=a6-::-SQRT ( 2 . 0* ( P2-P6 ) / ( RH06*k6 ) ) 

THE FOLLOWING EQUATIONS IN IMPLICIT FORM WERE INSERTED IN 
AN EFFORT TO CLEAR THE ALGEBRAIC LOOP PROBLEM. IT DID NOT 
WORK. 

THE NEXT LOGICAL STEP WOULD BE TO COMPACT ALL EQUATIONS 
SUCH THAT EACH VARIABLE IS EQUATED ONLY TO ITSELF OR A 
VARIABLE WHICH IS DIRECTLY COUPLED TO AN INTEGRAL. 

Q7=IMPL(Q7Z,0.001,Q7CC) 

Q7C=( (p5-P7-CQ57S«RH057*ABS(057(-0.5»RH07*K7Q7Q7/Q5s) . . . 

/ ( CQ57*MU57*LE57 ) )*2. -05 

Q7CC=Q7C 

q8=impl(q8z,o.ooi,q8cc) 

q8c=( (p6-p8-cq68s*rho68*q68«abs(q68)-o.5*rho8*k8*q8/a6s) . . . 

/ ( CQ68*MU68ttLE68 ) ) *2 . -Q6 

q8cc=q8c 

Q9=IMPL(Q9Z,0.001,Q9CC) 

Q9C=( (P7-P9-CQ79S*RH079*Q79»ABS(Q79)-0.5«RH09*K9*Q9*Q9/A7S. . . 

- ' / ( CQ7 9*MU7 9-”-LE7 9 ) ) *2 . - Q7 

Q9CC=Q9C 

Q10=IMPL(Q10Z,0.001,Q10CC) 

Q10C=( ( P8-P10-CQ8lS*RH08l*Q8l*ABS(Q8l)-0 .5*RH010*K10*Q10*Q10. . . 

/a8s ) / ( CQ81s-MU81*LE81 ) ) -“-2 . -q8 

Q10CC=Q10C 

Q11=IMPL(Q11Z,0.001,Q11C) 

Q11C=(P9-P11-CQ11S*RH011*Q11*ABS (Qll) )/ ( COll*KUll*Lll ) 
Q12=IMPL(Q127 , 0 . 001 , Q12C ) 

Q12C= ( P10-P12-CQ125*RH012*Q12*ABS ( Q12 ) ) / (CQl2-::-MUl2-::-L12) 
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S1?c-(I5-^-cq57sSrh657^7-abS ( Q5 7)-o.S.HH07«7^7«i7/*5S). 

q57 ( /(C057*M057*le57) 

q5?cc=q57c 

q 68 c=^P 6 -p 8 - CQ68s*RH068^ABS ( Q°8 ) • 5--RH08*k8-»q8*q8/ A&S) 

^ / ( cq68»KU68*LE68 ) 

q68cc=q68c 

^C=(P7-P9-ci^^ 

Q ' / ( CQ7 9*MU7 9*LE7 9 ) 

Q79CC=Q79C 

Q8lC=^P8-P10-CQ81S^RHO81i>ABS(Q81)-0 . ^RHOIO^KIO-^IO^IO. . . 

^ /a8s ) / ( CQ8 l-“-MU8 l-"-LE8 1 ) 

q81CC=q81C 

qil=cil*(ph-pi 2 ) 

B1I=CFA/P1+BI1 

B2I=CFA/P2+BI2 

bri=cfa/pr+bir 

b5i=cfa/p5+bi5 

b6I=CFA/p6+BI6 

B7I=CFA/P7+BI7 

B8I=CFA/p8+BI8 

B91=CFA/P9+BI9 

BlOI=CFA/PlO+3llO 

BlH=CFA/Pll+BIll 

B12I=CFA/P12+BI12 

b57I=cfa/p57+bi57 

B68l=CFA/P68+Bl68 

B79I=CFA/P79+ BI79 

B8lI=CFA/P8l+Bl8l 



->***DENSITY VALUES *SH 



hor-8 
DH1=8 
HO 2=8 

HO 5=8 

H06=8 
HO 7 =8 



OlO=J 

011 = 



..256E-05+8. 

.256E-05+8. 

.258E-05+8. 

.258E-05+8. 

.258E-05+8. 

.25&E-05+8. 

.258E-05+8. 

.258E-05+8. 

6.256E-05+8 

:8.258e-o5+8 

:8.256e-o5+8 

8.256E-05+8 

:8.256e-05+8 

:8.258e-o5+8 

:8.258e-o5+8 



ij.2 3E- 0 5* P R-"-BR I 
1^2 3E- 0 5’"' P 1'” 311 

k2 3E- 0 5'" P 2-”*E 2 1 
Ll 2 3 E- o 5* :: * P 5-- ; - B 5 1 
ii.2 3E-0 5*P8 -”*Bo I 
423E-05«P7*E7I 
42 3E- 0 5* P 8-::-b8 I 
423E-05*P8*B9I 

. 423e-o5-"-piO'-b >101 

. 423E-05-"-Pl 1 ‘"Blll 
. 1 l2 3 E” 0 5'"' P 1 2-U-B12 X 
,423E-05-~*p57*»*b 57I 
,i< 23E-05 -"-p68-*BB8i 
.423E-05--P? 9--B79I 
. 423E-05*p81*b81I 



*** VISCOSITY EQUATIONS **•* 

mu57=8 . 198e-06*exp ( P57*”-B57I ) 
MU68=8 . 198e-06*exp ( p68*B68l ) 
MU79=8. 198E-06*::-EXP ( P79*B79I ) 
MU8l=8 . 190E-O6*EXP ( p8l*B8lI ) 
WU11=8 . 198E-06-X-EXP ( P11*B11I ) 
MU12=8 . 198E-06-::-EXP ( P12*B12I ) 
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